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Pole mass of the W boson at two-loop order 
in the pure MS scheme 
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Department of Physies, Northern Illinois University, DeKalh IL 60115, 
Fermi National Accelerator Laboratory, P.O. Box 500, Batavia IL 60510 

I provide a calculation at full two-loop order of the complex pole squared 
mass of the W boson in the Standard Model in the pure MS renormalization 
scheme, with Goldstone boson mass effects resummed. This approach is an 
alternative to earlier ones that use on-shell or hybrid renormalization schemes. 
The renormalization scale dependence of the real and imaginary parts of the 
resulting pole mass are studied. Both deviate by about ±4 MeV from their 
median values as the renormalization scale is varied from 50 GeV to 200 GeV, 
but the theory error is likely larger. A surprising feature of this scheme is 
that the 2-loop QGD correction has a larger scale-dependence, but a smaller 
magnitude, than the 2-loop non-QGD correction, unless the renormalization 
scale is chosen very far from the top-quark mass. 
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I. INTRODUCTION 

The discovery [l|) 0 of the 125 GeV Higgs boson h at the Large Hadron Gollider (LHG) 
has completed the minimal Standard Model of electroweak symmetry breaking. Since the 
LHG has also not discovered any superpartners or other new fundamental particles, it is 
now more motivated than ever to perform precision analyses of the masses and interactions 
of the known particles of the completed theory. This paper concerns the complex pole mass 
[sl-l^ of the W boson. 
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calculated at 2-loop order. 

There have already been many studies HHSTj] that calculate contributions to the physical 
W boson mass, including all 2-loop order contributions and some QCD-enhanced effects at 
3- and 4-loop order besides. (These are reviewed in refs. [35|, |36|, for example.) Indeed, the 
accuracy of the most advanced of these calculations exceeds that of the present paper when 
it comes to predicting the IT-boson mass in terms of other measured quantities. However, 
the existing calculations have been done in on-shell or hybrid MS/on-shell schemes, or use 


expansions in small squared mass ratios, as in the case of ref. 2^, [2^. In this paper, I 


will provide a calculation that does not employ mass ratio expansions and uses a “pure” 
MS scheme, which means that the complete set of input parameters consists of only the 
renormalized running MS quantities 


'V, g, g\ A, yt, gs 


( 1 . 2 ) 


at a given renormalization scale Q. Here, v{Q) is defined to be the minimum of the radia- 
tively corrected effective potential in Landau gauge, which is now known to full 2-loop order 


38l | with 3-loop contributions at leading order in and yt [3^ , with Goldstone boson mass 


contributions resummed |40|, 1^. This allows v to be traded for the Higgs squared mass 


parameter m‘^{Q). The normalizations of v,rri^ and A are such that the Higgs potential is 


V = + A($1’$)l 


(1.3) 


and (<l>) = u/\/2, with a canonically normalized Higgs doublet held <h. 

In principle, the input parameters should also include the other quark and lepton Yukawa 
couplings, but these make only a very small difference in the present paper, as discussed be¬ 
low. In the pure MS scheme approach, all of the complex pole masses and other observables, 
such as the Fermi decays constant, are outputs, to be computed in terms of the quantities 
in eq. fll.211 . In practice, global hts to data may be used to obtain the relationship. In this 
paper, the input parameters of eq. fll.2p are all understood to be in the full non-decoupled 
(6-quark) Standard Model theory. Note that if the renormalization scale Q is chosen be¬ 
tween M\y and M*, the largest logarithms encountered in calculations of the physical masses 
of IF, Z, h, t will be at most ln(M|/M^) 1.5. 

It has been argued that the experimental vector boson masses My,exp as measured at 
colliders are related to the complex pole mass quantities by, approximately ^|, [^, [^: 


/pf2 _ jy^2 I p2 


(1.4) 


Numerically, this amounts to Mw,exp ~ Mw + ‘27 MeV in the case of the IF boson, assuming 
the Standard Model prediction for the width. Here, My is the the real part of the complex 
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pole of the propagator, while corresponds to what is sometimes called the “on-shell” 

mass. In the following, I will refer to Mw rather than Mw,exp- The current experimental 


43l is Mw,exp = 80.385 ± 0.015 GeV. 


value 

At the present time, the pure MS scheme is not quite competitive in numerical accuracy 
with the on-shell or hybrid schemes for the IV-boson mass calculation (although it is for 
the Higgs boson mass, which has been obtained to 2-loop order with the leading 3-loop 
corrections 


44 


45j). However, as the technology for loop calculations improves, it is quite 


possible that this will change. As a matter of opinion, I hnd the modular approach of 
the pure MS scheme to be conceptually simpler, and it can be easily extended to include 
contributions from new particles beyond the Standard Model, and the methods used can 
even be applied to other vector bosons (such as a W) in different theories. In any case, 
there is hopefully some value in being able to compare different schemes for the Standard 
Model observables, given their importance. 


II. W BOSON COMPLEX POLE MASS AT 2-LOOP ORDER 


In this section, I describe the calculation of the IH-boson complex pole mass. The calcula¬ 
tion reported here is restricted to Landau gauge, because only in that gauge has the effective 
potential been evaluated to full 2-loop order with leading 3-loop corrections, and this is nec¬ 
essary to obtain the relationship between the Higgs vacuum expectation value (VEV) and 
the Lagrangian squared mass parameter, used implicitly in the calculation below. However, 
the complex pole mass is a physical observable. It is therefore independent of the gauge 
hxing parameters 1^, as well as renormalization group invariant. 

In order to obtain the H^-boson complex pole mass, one hrst obtains, in terms of bare 
parameters in the regulated theory in d = 4 — 2e dimensions, the transverse self-energy 
function 


n(^) = (2.1) 

This is obtained by constructing the hh-boson self-energy function from the sum 

of all 1-particle-irreducible 2-point Feynman diagrams, and then contracting with — 
/p^)/{d — 1), where p^ is the external momentum and s = —p^, using a metric with 
Euclidean or (—,-f,-f,-f) signature. Factors of l/(167r^)^ are used to signify the loop order 1. 
Rather than including counterterm diagrams separately, it is more convenient and efficient 
to do the calculation in terms of the bare quantities: the VEV vb and the bare Higgs squared 
mass parameter m^, and the couplings Qb-i and A^, HtB, 933 , and then rewrite the results 
in terms of the MS quantities. 

The finite, renormalization-group invariant, and gauge-fixing invariant complex pole 
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squared mass can be written at 2-loop order: 


oZ. = »'b + + 7^^ [n< 2 '(tn) + n<‘i'(W'B)ni')(W'B)], (2.2) 




(Ibw 


where Wb = QbV^^/A. The bare quantities are then eliminated in favor of the MS renormal¬ 
ized parameters using: 
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(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 


to obtain in terms of the renormalized parameters. Here /i is the dimensional regular¬ 
ization scale. The MS renormalization scale Q is related to it by 


Q‘^ = 4:716 


( 2 . 10 ) 


where je is the Euler-Mascheroni constant. The conn 
exactly the same conventions as in this paper, in ref. 441 


erterm coefficients were listed, in 
, except for: 


„9 — 

‘'2,1 ~ 


‘' 2,2 


f/ + 3 / 9 = + - laV;, 

361 . 
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( 2 . 11 ) 

( 2 . 12 ) 


All of these counterterm coefficients can be obtained from the 2-loop beta functions and 
scalar anomalous dimension found in refs. BQ, j^; see for example the discussion sur¬ 
rounding eqs. (4.5)-(4.14) of ref. 


The procedure for the rest of the calculation is quite similar to that in ref. ^ , to which 


the reader is therefore referred for some more details, in a (perhaps futUe) attempt to avoid 
triggering the arXiv’s self-plagiarism detector. The Tarasov algorithm Mj is used to reduce 
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the 2-loop integrals to a basis set. The program TARCER that is often used for this 
purpose was apparently unable to handle a few of the necessary reductions in a hnite time, 
so I wrote a new Mathematica program RedTint implementing the Tarasov algorithm. (This 
program will be publicly released soon.) After expansion in e = (4 — d)/2, the Tarasov basis 
integrals were then written in terms of a set of basis integrals dehned and described in detail 
in refs. 


53 


54j |. The 1-loop basis integrals are: 


A{x), B(x,y) 


(2.13) 


and the 2-loop basis integral list is 

I{x,y,z), S{x,y,z), T{x,y,z), T(0,x,|/), U{x,y,z,u), M{x,y, z,u,v). (2.14) 

The arguments x,y,... are squared masses, and B, S, T, T, U, M also each have an implicit 
dependence on the external momentum invariant s = —p^, while A, B, I, S,T,T,U have 
an implicit dependence on the renormalization scale Q. The computer program TSIL 541 
can then be used for the efficient numerical evaluation of these basis integrals. TSIL uses 
Runge-Kutta integration of differential equations similar to that suggested in ref. 55|, and 
also includes relevant analytical results found in refs, j^. 


56 


After writing bare quantities in terms of MS quantities and expanding in e, the tree-level 
squared-mass arguments of the basis integrals used in the hnal result are: 


W 

z 

t 

h 


{g^ + g>V^, 

yy/2, 

2\v‘^ 


(2.15) 

(2.16) 

(2.17) 

(2.18) 


and 0 for photons and gluons. As in the Goldstone boson squared masses are eliminated 
by using the condition for the minimization of the effective potential after resummation. 


A- \v^ 


Ibvr^ 

^±g 


,'2 


yA{t) - 3XA{h) - y [3A(1T) + 2W] 
[3A(Z) + 2Z]} + ..., 


(2.19) 


as explained in section 4 of ref. ^ (see also ^ and 63|]). The same relation is used to 
eliminate from the tree-level Higgs boson squared mass, which appears as h rather than 
H = m? + 3\v‘^. In a future 3-loop calculation of the W (or Z) pole mass, the 2-loop version 
of eq. fl2.19p should be used; this can be found in eqs. (4.18)-(4.20) of ref. j40 |. 
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The 2-loop W boson squared pole mass is thus obtained, after hnally taking e —)■ 0, as: 


~ i^wMw 


= w + 


167r2^^ 


+ 


(IGtt' 
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(2),QCD 

W 


-j- /\ pi ,non—QCD 




( 2 . 20 ) 


where the right-hand side is a function of v, g, g', \,yt, g 3 ,Q, with all propagator masses 
expressed as W, Z, h,t, or 0. The list of 1-loop basis integrals used is 


= {A{h), A{t), A{W), A{Z), 5(0,0), 5(0, h), 5(0, t), 

5(0,Z), B{h,W), B{t,Z), 5(fT,Z)}, (2.21) 


while the list of necessary 2-loop basis integrals is: 


= {/(0,0,h), /(0,0,f), /(0,0,fT), /(0,0,Z), I{0,h,W), I{0,h,Z), 

1(0, W,Z), I{h,h,h), I{h,W,W), I{h,Z,Z), I{tA,Z), I{W,W,Z), 

S{h,h,W), S{h,W,Z), S{tA,W), S{W,Z,Z), T(h,0,0), T(h,0,f), 

T(h,0,fT), T{h,W,Z), T(f,0,0), T{t,0,h), T{t,0,Z), T(fT,0,0), 

T(Z,0,0), T{Z,0,t), T{Z,0,W), T{Z,h,W), T{0,h,W), Ti0,W,Z), 

5(0,f, 0, IT), U{0,t,h,t), f^(0,t,t, ^), U{h,W,0,0), U{h,W,0,t), 

U{h,W,h,W), U{h,W,W,Z), 5(fT,0,f,t), U{W,h,h,h), U{W,h,t,t), 
U{W,h,W,W), U{W,h,Z,Z), 5(fT,Z,0,0), U{W,Z,h,Z), U{W,Z,tA), 
U{W,Z,W,W), 5(Z,fT,0,0), U{Z,W,0,t), U{Z,W,h,W), U{Z,W,W,Z), 
M(0,0,0,0,0), M(0,0,0,0,Z), M(0,0,0,fT,0), M(0, 0, f, f, 0), M(0, 0, f, f, Z), 
M(0, 0, t, W, 0), M(0, t, W, 0, t), M(0, W, 0, Z, 0), M(0, W, t, h, t), 

M(0, W, t, Z, t), M(0, W, W, 0, W), M(0, W, W, h, W), M(0, W, W, Z, W), 

M(0, Z, t, W, 0), M(h, h, W, W, h), M(h, W, W, h, W), M{h, W, W, Z, W), 

M{h, Z, W, W, Z), M(W, W, Z, Z, h), M(W, Z, Z, W,W)}. (2.22) 


In each of the 5, S', T, T, U, and M integrals, the external momentum invariant is the 
tree-level squared mass, s = W. 

The 1-loop contribution to the pole mass is: 


AiJi = g^mV,,\g{b,t^V)+[NAnQ-\V,,\^-)+n^]f(0,0,W) 

2 ^ 'A(Z) 
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(2.23) 


where 


Nc = uq = riL = 3 


(2.24) 


are the numbers of colors, quark doublets, and lepton doublets in the Standard Model, 
respectively, and the fermion-loop function is 


f{x,y,s) = ■^\[{x-yf + s{x+ y)-2s‘^]B{x,y) + {x-y-2s)A{x) 


+{y — X — 2s 


+ {s — 3x — Sy)/9, 


(2.25) 


and the bottom quark mass and |Vtfep dependence has been included. The lighter quark and 
lepton masses can also be restored in the obvious way, by changing the 0 arguments of the 
function / in eq. fl2.23p and introducing additional Cabibbo-Kobayashi-Maskawa (CKM) 
mixing factors. Fortunately, however, the difference made by non-zero masses of 6, r, c,... 
and the presence of CKM mixing (assuming CKM unitarity and Vtb = 0.99914 ^) is less 
than about 1 MeV in both and for 50 GeV < Q < 200 GeV, and is much less for 
Q in the middle of that range, so those effects will be neglected for simplicity below. 

Note that 1-loop contributions involving 5(0,0), 5(0, Z) and 5(0, h) cancel, when the 
0 arguments correspond to Goldstone bosons and unphysical modes of the vector bosons in 
Landau gauge. This and similar cancellations in the 2-loop order part (mentioned below) 
are useful checks, as non-cancellation of such terms would have implied imaginary parts of 
the complex pole squared mass that do not correspond to any real decay mode of the W 
boson. 

The 2-loop QGD contribution is also simple enough to be written on a few lines in terms 
of the basis functions: 


A 


(2),QCD 

W 


= 939 


2 2 ( c 


24 


-4(t - Wf{2 + t/hF)M(0, 0, t, f, 0) 


+8(t - 21F)(1 + t/W)T{t, 0, 0) - (lOf + 81F)5(0, tf 

-{m/W + 56 + im/t)A{t)B{9, t) + (30tV^ + 42f - 12hF)5(0, t) 

-(40/hF + 24/t)A{tf + {m/W + 84)A(f) - 39hF + 17f/2 


-{hq - 1 ) 1 K {31 + 125 ( 0 , 0 ) + 81 KM ( 0 , 0 , 0 , 0 , 0 )} 


( 2 . 26 ) 















The remaining, non-QCD, 2-loop contributions, are much more complicated, involving a 
large number of terms. The form of the result is^ 


A 


(2),non—QCD 
W 




j<k 


(1,1) Al) ril) 

k -'7 -'fc 


+ (2.27) 


The coefficients cf''^ and and and are given in electronic form in an ancillary 
hie coefficients.txt provided with the arXiv source for this article. These coefficients 
are written exclusively in terms of the quantities W, Z, t, h, [by using eqs. (I2.15p - (l2.18p to 
eliminate g,g',yt, and A], as well as the hxed parameters Nc, uq, and ul. The latter can 
each be set equal to 3 in the Standard Model, but are kept general for checking purposes, 
and to tag the fermion loop contributions. 

It should be noted that the coefficients in the expression of the pole mass in terms of the 
basis integrals are not unique. This is because different basis integrals are related by special 
identities that hold when the s quar ed mass arguments are not generic. These identities 
include eqs. (A.15)-(A.21) of ref. 44|, and eqs. (A.14), (A.15), and (A.17)-(A.20) in ref. 64|. 

When setting s —)■ IT in eq. fl2.20p . one encounters singular behavior in individual terms, 
associated with photon lines attached to a IT boson propagator. In general, such potentially 
singular terms should cancel in the complex pole mass 2^. They are dealt with here by 
using expansions such as^ 


B{0,W) = I-A{W)/W+ {s-W)[l +A{W)/W-\n{W - s)]/W 

+(s - Wf[-l - A{W)/W + ffi(IT - s)]/W‘^ + 0{s - Wf (2.28) 


with ln(a;) = hi{x/Q‘^). Similar expansions of 2- loop basis functions that have thresholds or 
pseudo-thresholds at s = IT are carried out using the differential equations listed in section 
IV of ref. 53|, using methods similar to those found in After doing so, all pole and 
logarithmic singularities in s — IT that are found in individual Feynman diagrams cancel in 
the total eq. (I2.2np . an important check. 

Several other helpful checks were performed on the calculation. First, single and double 
poles in e cancel in This cancellation relies on agreement between the counter-terms 

(for X = v,g, g', A, yt, gs) as extracted from the (3 functions and Higgs scalar anomalous 
dimension in the literature, and the coefficients of divergent parts of the loop integrations 
performed here. Second, I checked that logarithms of G = m? + A^u^ cancel. This is required 
for the absence of spurious imaginary parts that could occur when the renormalization scale 
is chosen so that G < 0, and spurious divergences that could occur for G = 0. Third, I 


t Of the 78 coefficients for products of 1-loop integrals, 42 vanish. 

^ Eq. (12.281) is used to eliminate B{Q,W) everywhere, explaining its absence in eqs. (12.211) and (12.231) . 
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checked the absence of spurious imaginary parts of note that must be identically 0 
in the case ng = 1, = 0, because in the Standard Model the W boson can only decay to 

lighter fermion doublets. This checks cancellations between diagrams with Goldstone boson 
propagators and the corresponding Landau gauge vector propagator parts with poles at 0 
squared mass. Fourth, I checked that in each of the formal^ limits that the quantities W, Z, 
t, h, AW — h, AZ — h, At — Z, t — W, t + W — Z, or t + W — h vanish, the whole expression 
for is hnite and well-behaved, even though many of the individual 2 -loop coefficients 
in eqs. fl2.23p . fl2.26p . and fl2.27p are singular in one or more of those limits. This again 
reflects non-trivial relations between different basis integrals when squared mass arguments 
are not generic. Finally, the result for was analytically checked to be renormalization 
group invariant through terms of 2-loop order. In principle, this should be equivalent to the 
check of cancellation of 1 /e poles, but in practice it tests many intermediate steps of the 
calculation. This check is written as: 


0 = Q 


dQ 


= 

^pole 


d d 




X 


A 


.W 

’pole’ 


(2.29) 


where X = {g, g', A, yt, gs}, and 7 ,^ is the anomalous dimension of the Higgs field. It uses the 
derivatives of basis integrals with respect to the implicit argument Q, given in eqs. (4.7)- 
(4.13) of ref. and derivatives of the 1-loop basis integrals with respect to squared mass 
arguments, given for example in eqs. (A.5) and (A. 6 ) of ref. 14^. It also uses the beta 


functions and scalar anomalous dimension given in refs. 



A corresponding 


numerical check of renormalization scale invariance is performed in the next section. 


III. NUMERICAL RESULTS 


The numerical computation of s^jg given by eqs. fl2.20p - fl2.27p is accomplished using 
the program TSIL This requires only 13 calls of the function TSIL_Evaluate (which 
uses Runge-Kutta solution of coupled differential equations to obtain multiple basis integral 
functions simultaneously) as well as relatively fast evaluations of the integrals for which 
analytic formulas in terms of polylogarithms are known and incorporated in TSIL. 

For purposes of illustration, consider a benchmark set of input data: 


v{Mt) 

= 246.647 GeV, 

(3.1) 

9{Mt) 

= 0.647550, 

(3.2) 

g'm 

= 0.358521, 

(3.3) 


§ None of these limits are close to being realized in the real world. 
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X{Mt) 

= 0.12597, 

(3.4) 

yt{Mt) 

= 0.93690, 

(3.5) 

gsiMt) 

= 1.1666, 

(3.6) 


where Q = Mt = 173.34 GeV is the input renormalization scale. The top Yukawa coupling 
and strong coupling constant were taken from ref. 6^ version 4, and the electroweak gauge 
couplings were taken from ref. 36|]. The VEV v{Mt), which should minimize the radiatively 
corrected effective potential in the scheme used here, has also been chosen to approximately 
reproduce the experimental value of the Z boson physical mass, using a separate calculation 
(similar to the present one, and also in the pure MS scheme) that I plan to report on soon. 
The Higgs self-coupling A was simultaneously chosen so as to also obtain a Higgs pole mass 
of Mh = 125.09 GeV, using the calculation of 4^ as implemented in the program SMH 45|, 
at an optimal renormalization scale Q = 160 GeV. However, in the absence of a true global 
£t to available data, it is important to emphasize that the benchmark parameters chosen 
here should be viewed as illustrative, rather than as a prediction of Mw- 

The results for the renormalization scale dependences of and Tw obtained from 

in various approximations, are shown in Figures ITTI and To 

make the graphs, the input parameters v,g,g', X,yt,g 3 are run, using 3-loop beta functions 


66 l. l67l |. from the input scale Mt to the scale Q on the horizontal axis, and is re¬ 
computed at that scale. In the idealized case, Mw and Vw would be independent of Q if 
computed to sufficiently high order in perturbation theory. 

In Figure [3Tl the (green) dotted line is the tree-level result IV, which shows a severe scale 
dependence, due to the running of g and v. This is still large, but reduced, in the 1-loop 
result, given by the (red) short-dashed line. The majority of the remaining scale dependence 
is eliminated by including the QGD part of the 2-loop result from eq. fl2.26p as shown in 
the (blue) long-dashed line. The (black) solid line shows the full 2-loop result. Note that 
despite the large scale dependence of the 2-loop QGD correction, it is actually smaller than 
the 2-loop non-QGD correction in magnitude except for Q < 85 GeV, where the effect of 
ln(f) starts to become large. The 2-loop non-QGD correction is of order 40 MeV, but is seen 
to have a quite mild scale dependence. 

In Figure 13^ the (red) short-dashed line shows the running of Tiy computed at 1-loop 
order. Adding in the 2-loop QGD contribution, as shown by the (blue) long-dashed line, is 
a signihcant effect, but does not eliminate the scale dependence, which is mostly due to the 
electroweak 1-loop renormalization group running of g and v. However, including the 2-loop 
non-QGD corrections to greatly ameliorates the scale dependence, as it captures and 
compensates for most of the effect of running of g and v. 

For the range 50 GeV < Q < 200 GeV, the deviations of Mw and T w from their median 
values are both about ±4 MeV. For Mw, this is shown in close-up as the solid line in Figure 
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FIG. 3.1: The mass My/ of the W boson, obtained from the complex pole squared mass 
■^pole as a function of the renormalization scale Q at which is computed, 

in various approximations. The (green) dotted line is the tree-level result W, the (red) short- 
dashed line is the 1-loop result, the (blue) long-dashed line is the result from the 1-loop 
and 2-loop QCD contribution, and the (black) solid line is the full 2-loop order result. The 
input parameters v, g, g', X,yt, gs are obtained at the scale Q by 3-loop renormalization group 
running, starting from eqs. (IHTD-drel). Note that the usual Breit-Wigner mass Mu/,exp is 
about 27 MeV larger than Mw- 
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At Q= 173.34 GeV: 

V = 246.647 GeV, g = 0.647550, 
g' = 0.358521, ?i = 0.12597, 
y, = 0.93690, 93 = 1.1666 
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FIG. 3.2: The width Tpu of the W boson, obtained from the complex pole squared mass 
s|^je = ~ as in Figure [3Tj The (red) short-dashed line is the 1-loop result, the 

(blue) long-dashed line is the result from the 1-loop and 2-loop QCD contribution, and the 
(black) solid line is the full 2-loop order result. 
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FIG. 3.3: Close-up of the scale dependence of the mass M\y of the W boson, obtained from 
the complex pole squared mass — iTwMw, as in Figure [3d] The solid line is the 

full 2-loop order result, while the dashed line is the same, but after expanding the MS mass t 
(in the 1-loop part only) about T = (173.34 GeV)^ to first order, using eq. (13.7p . 

13.31 While this gives some lower bound on the remaining theory error (not counting the 
parametric errors in the inputs n, g, g', A, yt, gs), it is always questionable to assume a direct 
relationship between scale dependence and theory error. For another handle on the theory 
error, consider the following exercise. In the top/bottom 1-loop contribution, the running 
top mass t is used in propagators in the pure MS scheme. However, once the result has been 
obtained, one can expand t about any other value, for example the top-quark pole mass T. 
Doing so for the 1-loop contribution only is sensible, since t only appears in propagators, 
not vertex couplings, in the 1-loop order W boson self-energy. The relevant expansion is: 

/(O, f, W) = /(O, T, W) + {t- T) [A{T) -2W + {T + W)B{0, T)] jTW + 0{t - Tf. (3.7) 

If this expansion is extended to, say, 4th order in t — T, then the results are easily checked 
to be nearly indistinguishable from the original /(0,t, IF) without expansion. However, 
terminating the expansion at linear order in t — T, as in eq. fl3.7p . can be considered an 
alternative consistent 2-loop order result, if t — T is treated as formally of 1-loop order. 
This version of Mw is shown as the dashed line in Figure 1331 It clearly has a worse scale 
dependence, particularly at larger Q, where T — t becomes large. This suggests that the ±4 
MeV scale dependence of the original (solid line) pure MS calculation may be at least partly 
a fortunate accident. The two curves agree near Q = 77 GeV, where the running top-quark 
mass t equals the physical mass T. 
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IV. OUTLOOK 


In this paper I have reported the results for the complex pole mass of the W boson in the 
Standard Model in the pure MS scheme, with the vacuum expectation value, defined as the 
minimum of the Landau gauge effective potential, taken as one of the input parameters. The 
organization of input and output parameters is quite different from previous works that use 
the on-shell scheme or hybrid MS/on-shell schemes. The state-of-the art computations in 
these schemes, see respectively e.g. 31| and and references therein, probably both attain 
a better theory error than the pure MS scheme, for now. Moreover, a direct comparison of 
numerical results will need at least the corresponding results for the Z boson, which I hope 
to report on soon. Both results will then be incorporated into a publicly available computer 
code together with the Higgs boson mass code from 4^, l45 1. 


Refs. 


25 


26| and the very recent ref. 37| (which appeared as the present paper was being 


hnished) also used the pure MS scheme to compute the complex pole mass of the W boson. 
However, attempts at direct comparison are complicated^ by the fact that these papers used 
a different definition of the VEV, namely = —rn?/ A, rather than v that minimizes the 
nil radiatively corrected effective potential as made here (and, for example, refs. and^- 


36l|). The choice of using Utree requires including non-trivial tadpole diagrams, unlike the 


choice of expanding around v where the sum of Higgs tadpole diagrams (including the tree- 
level tadpole) simply vanishes. This means that already at 1-loop order, the expressions 
appear different. Compared to A|^^/(167r^) in eq. fl2.23p of the present paper, the sum of 
the bosonic contributions in eq. (B.2) of ref. 25| and the fermionic contributions in (B.2) of 
ref. differs by: 


9 


IQ'n'^h 


-2NM{t) + ^hA{h) + 3WA(W) + 2W^ + ^ZA{Z) + Z^ 


(4.1) 


This is simply because the tree-level terms are also different, namely /A in the present 


paper and g^vl^^jA in refs. l2y, l37|. To 1-loop order accuracy, the two expressions for the 


pole mass can easily be checked to be the same, by using eq. fl2.19p above, but establishing 
the connection at 2-loop order would require a somewhat non-trivial re-expansion using the 
2 -loop relation between and 

Note that, in general, expanding around Utree rather than v has the effect of making 




Also, refs. [25|,l2^ use expansions in 1/4 —sin^ 9w and Z/h and Z/t (in the notation of the present paper), 
which further increases the difficulty in making a direct comparison. 


t However, Ref. 


uses Feynman gauge instead of Landau gauge, so the VEV referred to in that paper will 


also not be the same thing as v in the present paper. Note that using v requires choosing a gauge-fixing 
prescription; choosing Landau gauge has the advantage that the effective potential is much simpler. 
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the perturbative expansion parameter be —rather than the usual —for the terms 
leading in the top mass. This can be seen in the presence of the hrst term in eq. (EH; 
in contrast, there is no gHA{t)/h term in eq. fl2.23p . As mentioned above as one of the 
checks, at two-loop order there is also no behavior like t^/h? or t^/h (or any other pole 
singularity in h, or W or Z) in in eq. fl2.27|) . As another example, see the 


discussion surrounding eqs. (4.34)-(4.40) in ref. where the terms of order 


yi 


lOvr^A 


in the relation between Utree and v are explicitly identihed for loop orders i = 1,2,3 in the 
limit ^ X in the case g = g' = 0. Not surprisingly, expanding around the radiatively 
corrected VEV leads to faster convergence than expanding around the tree-level VEV, at 
least formally, although both expansions should converge given enough loop orders, since 

is still numerically small. 

167r"*A 


It would clearly be useful to include the 3-loop contributions to W and Z complex pole 
masses in the pure MS scheme, so that theory errors can be made unambiguously much 
smaller than all relevant experimental errors. Here it should be remarked that it is not 
at all obvious that the parametrically QCD-enhanced contributions at 3-loop order will be 
the largest, especially considering that this was not the case at 2-loop order. A possible 
scenario is that the QCD-enhanced contributions will have the largest renormalization scale 
dependence, but not the largest magnitude, since this is what happened at 2-loop order. It 
seems feasible to eventually include all 3-loop contributions to in the pure MS scheme, 
although to do so without using mass expansions or approximations may require developing 
new methods for treating 3-loop self-energy contributions. 
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